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We consider gauge theories from the free evolution point of view, in which initial data satisfying 
constraints of a theory are given. Because the constraints are compatible with the field equations 
they remain so. We study a model constrained Hamiltonian theory and identify a particular structure 
in the equations of motion which we call the standard gauge freedom. The pure gauge subsystem 
of this model theory is identified and the manner in which the gauge variables couple to the field 
equations is presented. We demonstrate that the set of gauge choices that can be coupled to the 
field equations to obtain a, properly defined, wave-like formulation is exactly the set of wave-like 
pure gauges. Consequently we analyze a parametrized family of formulations of general relativity. 
The generalization of the harmonic gauge formulation to a five parameter family of gauge conditions 
is obtained. 
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/. Introduction. Field theories often have wave-like, 
or hyperbolic, degrees of freedom contained somehow in 
a set of variables, some of which are constrained, and 
some of which, the gauge fields, are not determined by 
the theory [l| . Physical states are equivalence classes of 
solutions related by a change of gauge. Solutions to the 
theory can be understood through properties of the equa- 
tions of motion, which consist of a mixture of the gauge, 
constraint and physical quantities. Unraveling this struc- 
ture in general may be hopeless. But if the gauge is care- 
fully chosen, say by taking the harmonic gauge in general 
relativity Q , then the full set of equations of motion may 
be rendered strongly hyperbolic [1, H|. This condition 
guarantees the existence of a unique solution to the initial 
value problem that depends continuously on the initial 
data, at least locally in time. As highlighted in Q "Ide- 
ally, one would like to exhibit a kind of hyperbolic skele- 
ton of the Einstein equations and a complete character- 
ization of the freedom to fix the gauge from which all 
hyperbolic reductions should be derivable. Instead, there 
are at present various different methods available which 
have been invented to serve specific needs," this ad-hoc 
characterization is unsatisfactory. Equations of motion 
for the gauge choice can be obtained in the absence of 
any coupling to the theory, which begs the question - 
what is this skeleton? In other words, what are the set 
of pure gauges that can be coupled to the theory to form 
a hyperbolic formulation? Since the basic characteriza- 
tion of a set of partial differential equations can be made 
in the linear approximation, we may start by directing 
our efforts there. We thus begin to address these issues 
in section II for a model linear constrained Hamiltonian 
system. In section III we examine conditions under which 
a formulation of the Hamiltonian theory is strongly hy- 
perbolic. In section IV we apply our findings to general 
relativity (GR) with a five parameter family of gauge 
conditions and obtain the generalization of the harmonic 
formulation to this family. 



II. A model theory with gauge freedom. Consider 
the equations of motion for the Hamiltonian density, 



H = 



IfdiqYfV*' Ft* \ f d jq 
2 \ p J \ F j M" 1 J \ p 

+ glCuV^didjq + glCM'M^dip , (1) 

with canonical positions and momenta (q,p). Every ma- 
trix is constant, M^ 1 is invertible and F % = f3 l I for some 
shift vector f3 l , with I the appropriate identity. Such a 
Hamiltonian can be obtained from that of GR by lineariz- 
ing Q and discarding lower derivatives. Variation with 
respect to the gauge fields (g q ,g p ) reveals the constraints 

n = C n V ij did jq = , M = Cm'M ] <>,l> = , 

which we will take to be first class and call the Hamilto- 
nian and momentum constraints respectively. 

Gauge invariance: The constraints generate the 
gauge transformation, 

q^q-=q-M- 1 C M U diilJ, 

p^p = p-V ij C^d i d j e ) (2) 

with unspecified fields and ip. We require that the field 
strength V^didjq and curl e I M~ 1 dip, denned by some 
square anti-hermitian matrices e l , are invariant under 
this transformation. Gauge invariance thus gives, 

{A H m) {i Cm 3) = C n V ij , V^M-'CW^ = , 
e ( i M- 1 V j VC n i = 0, 



for some matrices (A^vi) 1 , where the index parentheses 
denote symmetrization. Gauge invariance of the evolu- 
tion equations also implies that 



d t 6 = ?di6 + g q , 

d t i> = (A nM y%e + mil J + g p , 



(3) 



where g q , g p denote the change under and where here, 
and in what follows, /3 l is taken to include the appropriate 
identity. 
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Electric and magnetic degrees of freedom: Without 
loss of generality the matrix Cu has linearly indepen- 
dent rows, so we can decompose the potential matrix 
according to 



(4) 



with Cu = [CuC-h^]~ 1 ^ 2 C-h, an d the projection opera- 
tor _L-p defined implicitly by Q . In the absence of Hamil- 
tonian constraints we assume that 



Vp , and always 



that V$ 



t< (A VB )ei + C m ^ 1 {A V m)Cm j) for some 
Hcrmitian matrices {Ay b ) and (Avm ) ■ The electric and 
magnetic degrees of freedom are, 



E = V$a i d j q, 



B = e'M^dip. 



which, up to coupling to the constraints, form a closed 
subsystem and are gauge invariant. Such fields can be 
similarly defined in the absence of Hamiltonian con- 
straints. They are not used in the analysis that follows. 

Closure of the pure gauge subsystem: We call an equa- 
tion of motion for the gauge fields a gauge choice. Here 
we consider only evolution conditions 

dtg q = {A gq g q yd igq + (A gqgp ) l d igp + {A gqP )p, 

dtg P = {A gp g q yd igq + (A gpgp yd igp + (A gpq )%q . (5) 

We assume that {A gqP ) = ACu and (A gpq ) % = BCm 1 + 
C 1 C-hM for some matrices A,B and C l , a restriction 
which can be dropped by altering our arguments slightly. 
Assume that we are given a solution to the theory. We 
have already seen that the field equations are invariant 
under the gauge transformation @ . The pure gauge sub- 
system ([3]) is closed by substituting the gauge difference 
from @ into ([5]), taking g q — s- g q and g p — > g p . 

Free evolution on the expanded phase space: We are 
free to modify the dynamics of the model theory away 
from the constraint satisfying hypcrsurfacc in phase 
space, provided that the constraint subsystem remains 
closed. We define new constraints (O, Z) with the same 
length as {g q ,g p ) respectively. We couple the new con- 
straints to the gauge conditions ([3]) by parametrized ad- 
dition according to 

9t9g = {Ag q g q ) l d. igq + (Ag q g p yd igp + (A gqP )p 
+ (A 9q e)e, 

9t9p = {Ag^Jd^ + {Ag p g p yd igp + {Ag^f Drf 

+ (A gpZ )Z. 

Likewise for the equations of motion 

d t q = M- X P + F%q - M- x C M ^d igv + (A qe )e , 
dtp = \ j + F^ip - V'Cn <),<),,,,, + {A pZ ydiZ 

+ (A pn )H. 



We choose equations of motion for the new constraints 

d t O = fids + (AezfdiZ + (Aen)n , 
d t Z = {AzefdS + ftdiZ + {A ZM )M . 

The constraint subsystem is closed by 

d t U = (Ane^didjQ + p%n + (A nM )%M , 
d t M = (AMzfdidjZ + (A Mn )%H + p%M , 

with matrices 

(A H e) 13 =C n V ij (A q@ ), 

(A MZ ) ij = c M {l M-\A pZ y\ 

(A Mn y = C M % M-\A pn ). 

Natural choice of variables: The next assumption 
is that the variables can be appropriately broken up. 
For this we assume that for every unit spatial vec- 
tor s% the rows of Cu and Cm s = Cm 1s % are con- 
tained in the span of the union of the rows of V = 
CuV ss and W — Cm s M~ 1 , which each have them- 
selves independent rows, and furthermore that the con- 
tractions X — VCu^ and Y = WCm^ are invertiblc. 
With these conditions we can define 

Ce = —X~ 1 Cu , 

C 4 , = -Y- x Cm s + (Au M y s C g {M - C M ^Y- l C M s ] , 
L = I- V 1 [VV^V - W^WW^W , 

and the decomposition of d s q and p into gauge, con- 
straint, and physical degrees of freedom, 



d 2 s 9 = Cgp + (Aee)e , 9 S 2 V = C^d s q 
H = Vd s q, M = Wp, 

d s P q =.Ld s q, P p =±p, 



(A^z)Z, 



is invertible. The names here serve only to identify the 
relationship between the pure gauge and constraints. 

Principal symbol of a formulation: Once the gauge 
and constraint addition parameters are fixed we say that 
we have a formulation of the theory The principal sym- 
bol of a formulation in the direction is 



P s 



ps ps 



Pr 






P" 



(G) 



We assume that the constraint addition parameters are 
annihilated by the projection operator _L. This restric- 
tion can also be relaxed. The pure gauge sub-block, 



ps _ 



[3 s 

(Aum)^ s 
-(A gqP )V^ 








/ 





(Ag q g q Y 

-(A gpq yw^ {A gp9q y 



o 

(Ag q g p ) 
(Ag p g p ) S ) 
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is exactly the principal symbol of the pure gauge subsys- 
tem described after equation . The off-diagonal block, 



gc 



( o 

(A/,e) 



(A ez ) 



(A z ) 



(Am) 






\ 

{A^m) 


o / 



with sub-matrices, 



(Aez) = (A ge )(Aez) + C (A pZ ) s , 
(Aen) = (A ee )(Aen) - A" 1 + C e (A pn ) , 
(A^ e ) = (A 4 ,z)(A^ e ) - (A n M) u (A e) + C^(A 9qe ) , 
(A^ M ) = (A^z)(A ZM ) - y 1 - (Aum^CoCm^Y- 1 , 
(A e ) = (A 9qV )V\A m ) + (A 9q& ) , 

(A z ) = (A gpq yw\A ilZ ) + (A gv z), 

parametrizes the coupling of the gauge fields to the con- 
straints. The constraint violating sub-block, 



ps _ 



( P s (AezY 

(AzeY P s 

(A H e) ss 

V (A M zY 



(Aen) ^ 
(A ZM ) 
P s (AumY 
(AmhY P s J 



is exactly the principal symbol of the constraint subsys- 
tem. Finally the physical sub-block, 



ps 

1 T, 



_L V s 



_L A/" 1 



contains neither constraint addition or gauge parameters. 

Strong hyperbolicity: A necessary condition for 
strongly hyperbolicity is that P s has real eigenvalues and 
a complete set of eigenvectors for every s 1 . Strong hyper- 
bolicity is equivalent to well-posedness, that is existence 
of a unique solution depending continuously on the given 
data, of the initial value problem @, 0, 0] ■ 
III. Basic properties of theories with the standard 
gauge freedom. Consider the theory of the previous 
section. Then: 

Lemma. No formulation is strongly hyperbolic if the 
physical sub-block is not. 

Proof. Obviously a necessary condition for diagonaliz- 
ability with real eigenvalues of ([5]), for any formulation, 
is that of □ 

Lemma. A necessary condition for strong hyperbolicity 
of a formulation is that the pure gauge and constraint 
violating subsystems are strongly hyperbolic. 

Proof. We need to show that if the matrix © is diago- 
nalizable with real eigenvalues then this property holds 
for the pure gauge and constraint violating sub-blocks. 
A diagonalizablc upper block triangular matrix has diag- 
onalizable blocks on the diagonal [J, App. A]. Moreover, 



the set of eigenvalues of the full matrix is the union of 
the eigenvalues of the diagonal blocks. The lemma fol- 
lows. □ 

Definition. // for every strongly hyperbolic pure gauge 
there exists a choice of constraint addition parameters so 
that the formulation is strongly hyperbolic, we say that 
the theory has the standard gauge freedom. 

Theorem. Given a theory with the standard gauge free- 
dom, a pure gauge can be used to form a strongly hyper- 
bolic formulation if and only if it is strongly hyperbolic. 

Proof. The result follows trivially from the Lemmas and 
the definition of the standard gauge freedom. □ 

IV. Application to GR. The ADM Hamiltonian |[ 
for vacuum GR is Hadm = —ct H + 2 /3 l M,, with Hamil- 
tonian and momentum constraints 



H = R- KijK 13 



Mi = D J Ka — DaK . 



Gauge freedom in the non-linear regime: We take the 
freedom to be to choose coordinates x M = (t, x l ) on space- 
time; qualitative features of the model carry over. The 
constraints are obviously spatially covariant. Given an 
additional upper case time coordinate T with normal vec- 
tor N a such that N a = W(n a + v a ), with Lorentz fac- 
tor W and spatial boost vector v l then 

{N) H = W 2 H - 2W 2 M V , 

_L -WM; = WA'U + 2W 3 M vVi - W 3 Hvi , 

where is the projection operator into slices of con- 
stant t, and subscript v denotes contraction with the ve- 
locity v l . The electric and magnetic parts of the Weyl 
tensor [9(, form a closed subsystem, up to coupling to 
the constraints, and from the point of view of the lower 
case observer the spatial part of the upper case electric 
and magnetic parts are 

L -WEij = (2W 2 - l)Ekj - 2W 2 E v(lVj) + W 2 E vvllJ 
+ 2W 2 e\ {l B l)k , 



1 XN) B .. = W * B .. _ W 2 e k ij E kv - W 2 e k m E jk , 



which shows that if the fields vanish in one foliation they 
vanish in every foliation. 

Linearized pure gauge subsystem: The linearized pure 
gauge subsystem is [id ]. 



d t 6 = U- ipiD'a + ftdiO, 
dt^ = V l + aD l 6 - 6D l a + Cpi/)* 



(7) 



where 6 = -n a A[x% ip* = - _L' a A[x a ], U = A[a] and 
V = A[f3 1 ]. Under an infinitesimal change of gauge the 
perturbation to the metric and extrinsic curvature are 
given by the York equations ll| with a — > 8 and (3 l — > 
ip z , which can be used to close the linearized pure gauge 
subsystem once we act on the gauge condition with the 
perturbation operator A. 
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Free evolution in the expanded phase space: We ex- 
pand the phase space by constraints and Zj, and 
parametrize the equations of motion for the gauge by 

d t a = —g\c?K + g 2 ad l (3 i + ftdiCt + 2c 1 a 2 Q , 
d t J3 l = a 2 [<?37 fe V j + <? 4 7 Vl^-fc - gsad'a + fid,? 

+ 2a 2 c 2 Z l , (8) 

with gi > and g 3 = 2(173 + 34) > 0, and for the remain- 
ing variables by 

d t Kij = -DiDja + a[IUj - 2A',A„, + A,, A + £,A,, 
+ 2c 4 a9(jZ :( -) + i c 5 aj lj d k Z k + | c 6 a^ i3 H , 
<9 t 6 = c 7 aif + cga^Z 1 + £,36 , 
<9 t Zj = c 9 aMi + c w adiQ + CpZi . 

Strong hyperbolicity for non-linear and variable coeffi- 
cient systems is defined, with additional smoothness con- 
ditions, by linearizing and working in the high-frequency 
frozen coefficient approximation [7[. In this approxima- 
tion the Hamiltonian density has the structure of (fjl . 

Strong hyperbolicity of the pure gauge subsystem: The 
principal symbol of the linearized pure gauge subsys- 
tem Q with gauge choice ©, where one must ignore 
constraint addition, has eigenvalues ±yfgi , ±u±, with 

1v± = gi + g 3 - g 2 g5 

± V (gi +93- gig$) 2 - 4(51 - .92)33 • 

The subsystem is strongly hyperbolic if g 3 > and either 

i) - 7^ .92 < gi and g 2 gs < 9i - ^\JgT zr g 2 ~\fgl, + .93 , 

ii) . g 2 = and g 3 ^ gi or g 2 = 0, g 3 = gi and g 5 = 1 . 

the second clause of case ii) . is that of generically distinct 
eigenvalues colliding without loss of diagonalizability. 

Strong hyperbolicity of the constraint subsystem with 
vanishing gauge- constraint coupling: Choosing 

ci = gi , c 2 = g 3 , c 3 = c 5 = c 6 = , 

C4 = 2c 7 = c 8 = c 9 = l, cio = 2 ( 1 + — J , (9) 

guarantees both that the off-diagonal block of the prin- 
cipal symbol Pg C vanishes and that the constraint sub- 
system is strongly hyperbolic. The eigenvalues of the 
constraint violating sub-block Pq are ±^/cio, which are 
guaranteed to be real inside the class of gauges we are 
considering, and ±1 with multiplicity three. 

Strong hyperbolicity of physical sub-block: The phys- 
ical sub- block is diagonalizable with eigenvalues ±1, at 
least up to a trivial normalization. Assuming smoothness 
of the background implies the continuity requirement for 
strong hyperbolicity in every block. 



Discussion: The choice ^) is the natural extension 
of the harmonic gauge formulation Q to the family of 
gauge conditions If a gauge in which the contracted 
Christoffcl symbol is chosen to appear in the shift con- 
dition, i.e when 174 = —\g 3 , the constraint addition pa- 
rameters correspond to those of the principal part of the 
Z4 formulation [12j . Otherwise it differs in the constraint 
subsystem. 

V. Conclusion. Stimulated by in which the possi- 
bility of identifying every hyperbolic formulation of GR 
was suggested, we identified a particular structure in con- 
strained Hamiltonian equations of motion. We examined 
how pure gauge is inherited by a formulation of a the- 
ory. With this structure the set of strongly hyperbolic 
pure gauges are exactly those that can be used to form a 
strongly hyperbolic formulation, in-line with the expec- 
tation of the physicist. We expect that the results can be 
generalized to include elliptic gauges. It will furthermore 
be of interest to treat the initial boundary value prob- 
lem. We used our findings to investigate hyperbolicity of 
a family of formulations of GR, generalizing [l3[ to non 
Hamiltonian formulations. Open questions include those 
relating to long-term existence with different gauges. 
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